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Based on the previously proposed notions of action operators and of quantum integrability, 
frequency operators are introduced in a fully quantum-mechanical setting. They are conceptually 
useful because a new formulation can be given to unitary perturbation theory. When worked out 
for quantum spin systems, this variant is found to be formally equivalent to canonical perturbation 
theory applied to nearly-integrable systems consisting of classical spins. In particular, it becomes 
possible to locate the quantum-mechanical operator- valued equivalent of the frequency denominators 
\^ • which may cause divergence of the classical perturbation series. The results established here link 

| the concept of quantum-mechanical integrability to a technical question, namely the behaviour of 

. specific perturbation series. 
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^ ! I. INTRODUCTION 

i— i ! 

■ In classical mechanics, the distinction between integrable and nonintegrable Hamiltonian systems is clear-cut. On 
the one hand, for a system to be integrable, the existence of N smooth, involutive, and functionally independent 
constants of the motion is required. Then, its 2A-dimensional phase space is guaranteed to be foliated into N- 
dimensional tori, and the resulting almost-periodic motion of the system is very simple. On the other hand, a 

, nonintegrable system may explore a (2N — l)-dimensional region of phase space in the course of time if (except for 
the energy) there are no invariants which restrict the motion to lower-dimensional manifolds. Accordingly, the time 
\ evolution of a nonintegrable system is extremely complicated . 

Much is known about the properties of quantum systems obtained by quantizing classically integrable and noninte- 
grable systems, respectively M. Nevertheless, no rigorously established and generally accepted features have emerged 
which would allow one to unambiguously distinguish integrable and nonintegrable quantum systems in analogy to 
classical mechanics. While the concept of integrability is important in classical mechanics (for example, the KAM 
theorem deals with the influence of a perturbation on integrable systems), it appears, for the time being, to be much 
less fertile in quantum mechanics. 

A new concept of quantum integrability has been proposed in || which is based on quantum-mechanical action oper- 

■ ators without any reference to classical mechanics. Incorporating criticism of earlier concepts of quantum integrability 
|Q| , this notion relies only on the algebraic structure of a quantum system provided by the commutation relations of the 
basic operators. It is the objective of the present paper to explore the features of nearly quantum integrable systems: 
a perturbation is added to a system which qualifies as quantum integrable in the above sense, and a method is given 
to successively remove the disturbance. By using the quantum actions in a framework of Lie transforms it becomes 
possible to lay open strong formal analogies to the corresponding classical perturbational approach. In particular, with 
frequency operators being defined for integrable quantum systems, it will become obvious that quantum-mechanical 
"counterparts" of the classical frequency denominators show up in this perturbational treatment. The results of the 
present work are formal, i.e., questions of convergence have not yet been dealt with. Fortunately, the structural 
similarity of the results to classical mechanics allows one already to imagine possible scenarios. 

Two general remarks are in order before turning to the detailed presentation. The investigations are carried out for 
spin systems but the generalization of the results to particle systems is not expected to provide substantial difficulties. 
Then, it should be pointed out that no part of this work relies on semiclassical approximations. 

In brief, the paper is organized as follows. In Section 2, the notion of quantum integrability used here is recapitu- 
lated. Frequency operators are introduced in Section 3, and they are recognized to naturally generalize the classical 
frequencies known from integrable systems. The main part of the paper, Section 4, consists in formulating the theory 
of Lie transforms for nearly quantum integrable systems. Applications to systems with one and two degrees of freedom 
follow in Section 5. Finally, the results are summarized in Section 6, and connections to other work are discussed. 
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II. THE ALGEBRAIC VIEW OF QUANTUM INTEGRABILITY 



It will be useful to summarize the algebraic view of quantum integrability in order to motivate the investigations 
of the subsequent chapters. Only the essential points will be mentioned; details have been elaborated on in Ref. S. 

Let us outline the concept of integrability with a view on the algebraic structure of classical mechanics. For a system 
of N degrees of freedom there are canonical variables q = (jji, . . . , q?j) and p — (pi, . . . ,pjv)- They provide a basis 
for the algebra of (smooth) phase-space functions f(p, q) since each such function can be approximated arbitrarily 
well by summing products of Pm Qn 

with appropriate coefficients. They can be multiplied by real coefficients, and 
one can add them; any two phase-space functions define a third one via the Poisson bracket turning the set of all 
phase-space functions into a Lie algebra. This framework is the same for all classical mechanical systems. A particular 
physical system is described by selecting one specific phase-space function and by calling it the Hamiltonian. Its role 
is to generate the time evolution of the system by defining a flow in phase space. Integrability comes in as follows. 
Typically, the Hamiltonian function will depend on all components of momentum p and coordinate q. However, it is 
possible to perform a canonical change of basis, (p,q) — ► (p',q'), as a result of which the Hamiltonian will turn into 
a different function if expressed in terms of the primed variables. It may happen that such a change of basis leads to 
a particular form of the new Hamiltonian function, namely that it depends on only half of the number of variables, 
N. One has found a particularly convenient basis of the algebra because it becomes obvious that the flow in phase 
space decomposes into N decoupled flows. A system with such a trivial flow is called integrable. As is well-known, 
this situation is very rare: typically no decoupling basis exists, and the phase-space flow is extremely intricate. 

It is straightforward to rephrase this line of thought for quantum-mechanical systems. Now the basis of the algebra 
of operators is given by p = (p\, . . . ,Pn) an d Q — (<Zi, • • • j Qn)i an d each operator in the algebra can be represented by 
an appropriate sum of products of p ni q n - The commutator of two operators defines a third one, and Jacobi's identity 
holds in general: one is dealing with a Lie algebra of operators. If one particular operator of the algebra is chosen as 
Hamiltonian, a quantum system and its time evolution are determined. A generic Hamiltonian operator will depend 
on all components of p and q. In special cases, however, a unitary transformation exists which allows one to express 
the Hamiltonian as a function of only half of the number of basic operators, N . By analogy, the system will then 
be called quantum integrable. Part of the purpose of the present work is to show that, in an appropriate sense, the 
corresponding quantum-mechanical flow is trivial, and that this situation is not generic. 

It should be noted that in the context of quantum mechanics the term "algebraic" will be understood as "inde- 
pendent of the representation." The idea is to only use relations which have a counterpart in the classical algebra. 
For example, within each finite-dimensional matrix representation of the spin algebra there exist relations between 
different powers of spin matrices which are specific for the representation at hand. Such relations do not exist in the 
classical algebra and, thus, they should not be used. An instructive discussion of the common algebraic structure 
underlying classical and quantum mechanics has been given by Falk ||. 

For convenience, the definition of quantum integrability for spin systems proposed in J3| will be reproduced. The 
description of a quantum spin system with TV degrees of freedom is based on N copies of the one-spin commutation 
relations, 

[Sja.Sjip] = ihdjj* ^2 £ af) 1 Sj 1 , j,f = 1,2, . . . N. (1) 

Here is the criterion for integrability of quantum spin systems: 

QJ: A given TV-spin Hamiltonian if — if (Si, . . . , Sjv) is quantum integrable if there exists a unitary transformation 
W(Si, . . . , S^v) which converts the spin operators Sj, j = 1, . . . , N, into new spin operators, 

( S i =WS i Wt = l S i (S l! ...S JV ), i = l,...,iV, (2) 

such that the Hamiltonian turns into a function of N quantum actions Jj = Sj Z : 

H(S 1 ,...,S N ) = H(J 1 ,...,J N ). (3) 

The spectra of the operators jj, j = 1, . . . , N, in each (2s + l)-dimensional representation T s consist of uniformly 
spaced levels: 

+s 

Jj = ^2 \---,m j ,...;s)m j h(...,m j ,...;s\, j = l,...,N. (4) 

Ulj — — S 
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If the transformation (|3j) has been achieved, the eigenvalues of the Hamiltonian can be read off immediately, and its 
eigenfunctions are known as well. In spite of reference to the representation r s , Eq. is an algebraic statement 
since it is required to hold in all representations, s = 1/2, 1, 3/2, ... As will become clear in Section 4, condition (||) 
defines a program which can be carried out in strong analogy to classical mechanics. 



III. FREQUENCY OPERATORS 



Frequency operators, to be defined for quantum integrable systems in a straightforward way, are introduced in 
this section. The main objective of Ref. [g| has been to provide numerical and analytical support for the existence 
of quantum- mechanical action operators, and to use them as building blocks for quantum-mechanical integrability. 
These arguments will not be repeated here — instead an independent approach will be presented which illustrates the 
idea of action operators from a new perspective. To arrive at the definition of frequency operators will then be a small 
step. 

As a starter, a classical spin system with one degree of freedom is considered. Its Hamiltonian is assumed to depend 
on a single spin component only, S z , say: 



The complex quantities 



H = H{S Z ) . (5) 



S± = S x ± iS y (6) 



are classical counterparts of ladder operators for a quantum spin. In these variables the Poisson brackets for a single 
spin read 

{S Z) S±} = ±iS±, {S + ,S^} = 2iS z , (7) 
as obtained from the replacement [ , ] — > — ih{ , } in (Q). The equations of motion are given by 

§ = {#,S}, (8) 

and for the Hamiltonian ([s]) they read explicitly 

^- = {H,S z } = 0, (9) 

^± = {H,S + }=iS + u,(S z ), (10) 

the equation of motion for S~ being complex conjugate to that of S+. When expressing the spin as S = 



(y S 2 ~ V% costp, \J S 2 — p^shnp^p^) in terms of canonical coordinates with {p v ,tp} = 1, it becomes obvious that 

the third component, S z , plays the role of an action: it is both a canonical variable and a constant of the motion. 
The phase-space functions S± are proportional to the exponentiated angle variable. In view of later developments it 
is useful to stick with the periodic functions of time, S±(t), instead of using the canonical variable tp. The frequency 

= ^ (ID 

is completely determined by the Hamiltonian and thus depends on the action S z only. The time evolution of integer 
powers k = 2, 3, . . . of S+(t) is obtained from ( [To| ) as 

dS k 

' - {H,S%} = iS%u k (S s ), (12) 



dt 



and similarly for S-. The notation 



u k (S M ) = ku(S t ) , (13) 
used here will simplify comparison with quantum-mechanical expressions. 
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Consider now a quantum-mechanical spin system with one degree of freedom, i.e., one copy of the commutation 
relations (Q) only occurs. The Hamiltonian of the system is assumed to depend on the third component of the spin 
operator S only, 

H = H(S Z ) . (14) 
Again, ladder operators are useful: S± = S x ± iS y , leading to commutators similar to (0): 

[S„S±] = ±hS± , [S+, S-] = 2hS z . (15) 
Heisenberg's equations of motion for the quantum spin are 

? = i^.S]= (16) 



dt h l 



and for the Hamiltonian ( |l4|) they read 



dt h 

and the equation for SL = (<S+)^ is the adjoint of the last equation, 

dS- i 



' /,Sf = M,S+]=iS + u,(&), (18) 



[H, S-] = -iuj(S z ) S- . (19) 



dt n 

The selfadjoint operator w = co(S z ) depends on the z component of the spin operator S only; explicitly, one has 

lo{S z ) = i (H(S S + ft) - H(SS) , (20) 
providing thus the definition of a quantum-mechanical frequency operator. Its form is suggested by requiring that 



the operator S+ stand on the left of everything else in Eq. fllq). The derivation of (20) hinges on the fact that the 



Hamilton operator depends on the operator S z only because in ( |lq ) one can then use the relation 

[f(S z ),S%] = S\ (f{S z + kh) - /(&)) = Sl(cxp(khd z ) - l)f(S z ) , (21) 



where d z denotes the (formal) derivative with respect to the operator S z . An algebraic proof of (J21J) is given in 
Appendix A for functions f(x) which have a power-series expansion in x, and another one for arbitrary functions f(x) 
which holds in each representation. A different but equivalent definition of the frequency operator can be given if one 
moves S + to the right in Eq. ([l^). Either convention is unambiguous; in the following, the version occurring in JTs| ) 
will be used. 

Comparing the classical and the quantum-mechanical equations of motion, ( |l0| ) and (|l8|), it becomes plausible to 
call the operator S z an action operator: it is both a constant of the motion and (the closest analog of) a canonical 
momentum in the spin algebra. It is reasonable to consider the Hamiltonian ( |l4| ) as quantum integrable since it 
depends on nothing but an action operator. This also fits well with conceiving uj given in (EOT) as frequency operator 



because (f20j) is the discretized version of ( |ll[ ) 



On this basis, the notion of quantum-mechanical integrability as defined in QJ makes sense: a system with Hamilton 
operator TL will be called quantum integrable if a unitary transformation U can be found such that the Hamiltonian 
turns into a function of commuting z components of the different spins only while the commutation relations (|l|) 
remain invariant. As a consequence, the equations of motion in the Heisenberg picture strongly resemble those of 
classically integrable systems. In addition, frequency operators are defined in analogy to the classical frequencies both 
of which are determined completely by the Hamiltonian. It is not for the first time that operator-valued analogs of 
the classical frequencies do occur: as early as 1925, Dirac introduced a similar concept for particle systems ||. 

It is important to look at the time evolution of powers of S+ (and SL). One finds 
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dS k 



dt h 

where now 



-[H,S k + ]=iS k + oj k (S z ), (22) 



Since generally 



u Jk (S z ) = l(H(S z + kh)-H(S z )) . (23) 



£ fe + , (24) 



a difference to the classical relation (|13|) arises which will have interesting consequences for the perturbative approach 
to be studied later on. When formally expanding both sides of ( pi| ) in powers of h, however, the leading order reflects 
the classical relation. 

The generalization to systems consisting of N spins, Sj,j — 1, . . . , N, is staightforward. The components of each 
individual spin fulfill commutation relations (||), and operators associated with different spins Sj a , Sj'p commute. The 
Hamiltonian now will depend on all spins, H = -ff (Si, . . . , Sat) and the equations of motion read 

dS- i -~ 

—rjr- = -AH, S ja ], a = x,y,z. (25) 
For a quantum integrable iV-spin Hamiltonian there are N frequency operators generalizing the expressions ( p0| ) and 

Lo { k 3) (S z )^^H(S lz ,...,S ]Z + kh,... 1 S Nz )-H(S z )^ , (26) 
and the time evolution of the ladder operators is simply given by 



dS k 



-^ = -[H,Sf + ] = iSf + QV> , fc = l,2,..., (27) 
and their adjoints. 



IV. PERTURBATION THEORY 



Having established a concept of integrability for quantum systems, the program to be carried out is as follows. In 
classical mechanics, a small perturbation is added to an integrable system and it is investigated whether, by successive 
canonical transformations, the perturbation can be removed completely; this would render the system integrable. In 
consequence of the KAM theorem, however, integrability is destroyed by a generic perturbation: any perturbation 
series is doomed to diverge in finite fractions of phase space. A quantum-mechanical implementation of this approach 
aims at transforming a perturbed quantum integrable system into an integrable one by applying a sequence of unitary 
transformations. Lie transforms are a convenient tool here since they are adapted to the algebraic structure of the 
theory. Moreover, they facilitate comparison with canonical perturbation theory, the relevant formulae of which are 
gathered in Appendix B. 

A quantum system consisting of N spins with constant length S? is assumed to be described by a Hamiltonian 
operator of the form 

H £ =H{S,e) = H (S z ) + H(S,e), (28) 

where Hq is a quantum integrable system, depending only on the z components of the spins (if not indicated otherwise, 
S denotes the collection (Si, . . . , Sat) of operators). 

The perturbation H (S, s) may depend on all components of the N spins and it contains a real expansion parameter 
e such that the perturbation vanishes for e = 0: H(S, 0) = 0. Furthermore, it is assumed that the unperturbed 
Hamiltonian Hq is not degenerate. In the following, a possible dependence of operators on the invariants 5| will not 
be made explicit. 
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A. Unitary transformations 



To begin with, it is useful to study the effect of unitarily transforming the spin operators by an operator U £ — U(S,e) 
with U(S,e — 0) — 1. New spin operators S' (throughout this paper the prime does not denote a derivative) are 
introduced by 



= U e $) §[/](§), 
and the induced transformation of the Hamiltonian reads 

H(S,e) = H'(S' e ,e), 



(29) 



(30) 



that is, in terms of the primed operators, S^, the Hamiltonian is given by a new function H' . By using the inverse of 
(^) in the form 



relation (30) can be written as 



s = ut@)s' e u e @) = u2@ e )s' e u$' e ), 



U'J(S' £ )H(S' e ,e)U'M) = H'(S' £ ,e). 



(31) 



(32) 



This equation holds for any unitary transformation U e . In view of the criterion for quantum integrability, QJ, the 
new Hamiltonian "H — H' is required to depend on the z components of the primed spin operators S = S' only, so 
that (30) assumes the special form 



H(S,e)=H(S„e). 

Defining U £ (S) = U £ (S), Eq. © turns into 

Z4(«S) H(S,s)U e (S) = H(S„e) 



(33) 



(34) 



In order to distinguish quantities related to a general transformation as in (^) from those associated with (|33h, curly 
symbols such as U and Tt will be used. It will be convenient to rename the curly operators S by S in Eq. (p4[). 



B. Perturbation expansion 

The construction of the unitary transformation U £ achieving ([m]) is at stake now. Under three assumptions it can be 
reduced formally to the successive solution of an infinite set of hierarchical equations determined by the unperturbed 
Hamiltonian Hq(S z ) and the perturbation H(S, e). First, it is assumed that the perturbation in J2g| ) can be expanded 
in a power series of the parameter e: 

oo 

J?(S,e) = £y #„(§). (35) 

When writing the operator IA £ as a Lie transform, 

W(S,e) =exp(-TO(S,e)/7i), (36) 
a similar expansion is required to exist for the hermitean operator u £ : 

oo 

u(S, E ) = X> P %.(S). (37) 
P =i 

The operator u £ is called the generator (or generating operator) of the unitary transformation IA £ . Finally, it must be 
possible to expand the new Hamiltonian operator, Ti., in powers of s: 



G 



W(&,e) = Xy«p(& 

p=0 



(38) 



Plugging these expansions into requirement (34) and collecting terms multiplied by equal powers of the parameter e, 
one obtains the following set of equations: 






= H 


- Ho , 






h 


= H 1 


- Hi , 








= H 2 


- ^[ui,Hi + Hi] 






i ^ ^ 
j[H ,u p \ 


= Rp 


Ml,.- - . . 


■ j Hp, ri\, . 


■ ■ i Hp-i) — H. 



(39) 

For convenience, each equation has been solved for the commutators of the unperturbed Hamiltonian, Hq, with the 
generator of the highest index, u p . Eqs. (39) constitute a nested hierarchy: at the pth. level the quantities Tip and 



Up can be determined if all quantities with smaller indices have been found; the operator R p contains only quantities 
know from the beginning or determined when solving the equations of the lower levels. Such a structure is typical for 
perturbation theory formulated in terms of Lie transforms. It is possible to discuss the method to solve the equations 
( |39| ) in full generality For simplicity, the procedure is carried out here for a system with a single degree of 

freedom first; generalization to two (or more) spins is then straightforward. 



C. One degree of freedom 

For a single-spin system with Hamiltonian 

H(S,e) = H Q (S z )+H(S,s), (40) 

the unknowns u p and Hp are determined from the nested hierarchy ( |39| ) as follows. 

Solving the first equation of the hierarchy is particularly simple because it states that the old and the new Hamil- 
tonian, Hq and Ho, respectively, coincide. By assumption, Hq depends on the z component of S only, and, thus, Hq 
has this property, too, while S and S are identical to order e . This is consistent as is easily seen from the limit of 
vanishing perturbation, e = 0. 

The remaining equations are all of the same type: 

^[H (S z ),u p (S ± ,S z )} = R P (S ± ,S Z ) - H P (S Z ) . (41) 
The operator R p can be written as a power series in the operators S± (cf. Eq. ( |84| ) in Appendix C): 

oo 

R P (S ± ,S Z ) = R p ,o(S z ) + (slR+ k (S z ) + R- k (S z )S k _) , (42) 

with uniquely defined operators Rp k — (R~ k y and R p .q which depend on S z only. 
On the other hand, expanding u p (S±, S z ) as in 

oo 

u p (S ± ,S z ) = Up. (S z ) + + u p,k&)S-) > (43) 

fe=i 

one obtains for the commutator in (|41 
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l - [H , u p ] = ]T (S h +u+ k iu k + (^iQ k )u^ k S k l) , (44) 



k=l 



where equation ( |22| ) has been used. The order of the operators k and the frequency operator of the unperturbed 

system, Qk, is immaterial since they depend on S z only. Obviously, there is no term depending on S z alone; hence, 
for (f4l| ) to hold one must have 



Furthermore, the operators ut k follow from comparing (f42|) and (p4|): 



H P {S Z ) = Rp }Q {S z ) . (45) 

fv"\l 1 i~\nr fv/^TVi Cf 

u p (S ± ,S z ) = u pfi (S z ) + V feS^ + , (46) 

k~i\ luJk ~ luJk J 

where the inverse of the operator is well-defined only if the operator oj k does not have a zero eigenvalue; this, 
however, is guaranteed by the assumption that the unperturbed Hamiltonian has no degenerate eigenvalues. Therefore, 



an operator u p solving Eq. (41) has been found in terms of R p and u) k . Its 'diagonal' part u p0 (S z ) (corresponding to 
a phase shift when acting on eigenstates of S z ) remains undetermined — for simplicity, it will be chosen equal to zero 
for all p. As a result, the possibility to solve the hierarchy of equations (|39| ) has been established. The similarity of 
Eq. ( f46| ) with the classical formula, Eq. ([79]), is striking. 

It is important to note that the expression (|4^) for the generator u p exhibits factors inversely proportional to the 
frequency operators uik(S z ), as is familiar from classical perturbation theory (cf. Appendix B). As can be read off 
from the structure of the hierarchy (|39"|), the operators R p are linear in u p -\: thus, they contain terms proportional 
to l/tDjjp 1 . Solving ( |46| ) for the generating operator of pih order, u p , one finds it to be proportional to the inverse 
of the pth power of frequency operators. The structurally similar classical expansion is bound to converge to a well- 
defined canonical transformation since all classical systems with one degree of freedom are integrable. By analogy, it is 
expected that the indivdual generating operators u p will add up to a sensible expression for the unitary transformation 

IA S . If classical systems with two degrees of freedom are studied, the impact of a perturbation is known to be disastrous 
in most cases. 

D. Two degrees of freedom 

For a nearly-integrable quantum system with two degrees of freedom 

^(Si,S a ,e)=flo(Si*,ft.)+ff(Si,S 2 ,e) ) (47) 

the transformed Hamiltonian operator should depend on two new spin components only, S\ z and S2 Z - As before, the 
lowest order of the hierarchy ( |39| ) is easily solved by defining 7^o to be identical with the unperturbed Hamiltonian, 
Hq. For p > 1, the equations to be solved read: 

j[Ho{S\ z , 622), u p (Si±, S2±, S\ Zl S2 Z )\ — R P {Si±, S2±,S\ Z , S% z ) — H P {S\ Z , S^z) . (48) 

Using expansions for the operators R p and u p (cf. Appendix C) in analogy to ( pl2| ) and ([43|), respectively, one finds 
for the generator 

t; (a a\- ( ok oi ^M^ 2 ) 

Up{S±,a z ) - ^ \ °1+ J 2+ - (!) - ~ — (2f 



(fc,<?)^(0,0) 



i{u k L, {s lzl s 2z + en) + 4" (s lz ,s 2z )) 



+Sf + m J' L ^ S l 2 _ + h.c. . (49) 

Again, the generator u p is determined up to an arbitrary function of S 2 only which is chosen equal to zero. When 
expressing the denominators in terms of the unperturbed Hamiltonian, it becomes obvious that both spins enter 
symmetrically in (f49|): the first one reads 
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~ (h (s u + kh, s 2z + en) - h (s 1z J 2z )) , (50) 

and the second one is equal to 

A« = ~ (Ho(5i* + fca, S 2z ) - flo(Su, S 2z + thj) . (51) 

Eq. (^9|) belongs to the main results of the present paper. For systems with two (or more) degrees of freedom the 
generators contain terms proportional to the inverse of linear combinations of frequency operators. The expansion 
( [49| ) makes sense only if the denominators Aj^ do not have a zero eigenvalue, which will be assumed from now on. 
When dealing with the convergence of the perturbation series, it will be necessary to proceed in two steps: first, the 
existence of each generator u p has to be shown, and subsequently it has to be checked that the sum of all u p leads to 
a sensible result . 

V. EXAMPLES 

The pertubation theory as developed in the previous section is applied here to explicit examples. Systems with 
one and two degrees of freedom will be studied in order to obtain more detailed insight into the structure of the 
denominators. 

A. One degree of freedom 

The Hamiltonian to be studied is defined as follows 

= aS z + \0l + |(§ 2 - Si) + 101 + Si) , (52) 

the nonzero constants a and j3 being real numbers, and the quantum integrable part, H = aS z + f3S 2 /2, is assumed 
to have no degenerate eigenvalue. The perturbation H(S, e) — sS 2 /2 is particularly simple because it consists of just 
one single term being proportional to e. 

It is straightforward to calculate the frequency operators according to (^3|) : 

Q k = k(a + f3S z ) + ^(3k 2 h, fc = l,2,... (53) 

The first term on the right-hand-side defines a function ujk of the operator S z which coincides with the expression for 
the classical frequency: u)k{S z ) — kuj(S z ) = k(a + \(3S Z ). Due to the discretization of the derivative an additional 
term proportional to fi shows up in (m) . 

According to (|45| ) the new first-order Hamiltonian is equal to the diagonal part of the perturbation: 

H 1 (S z ) = l(S 2 -S 2 z ), (54) 

as can be read off from Eq. (|5^). When expanding the off-diagonal part of the perturbation the only nonvanishing 
coefficients are found to be 

leading to the following expression for the first-order generator 

ur = UsljL + ^sA , (56) 

8 V IU 2 -ILU2 J 



with LJ2 from (53). Since the energy eigenvalues of the Hamiltonian Hq in ( |52| ) have been assumed not to be degenerate, 
the expectation value of the denominators in Eq. ( [S6l) will not have zero eigenvalues. Consequently, there is no 
principal obstacle which would destroy the convergence of the perturbation series for the unitary transformation U, 
similar to classical mechanics. 
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B. Two degrees of freedom 



Consider a two-spin Hamiltonian ( f47| ) with 

H(St, S 2 , e) = A(S U ) + / 2 (^ 2z ) + 7^1,52, + S(Si, S 2> e) , (57) 



where 



/■,■(&) =a j x+h j x 2 , j = 1,2. (58) 



The quantum intcgrable Hamiltonian consisting of the first three terms on the right-hand-side of Eq. ( p7| ) is quadratic 
in the z components of the spin. With a view to the structure of the denominators present in the generating operator — 
which are completely determined by the frequencies of the unperturbed system as seen from Eq. ( |49| ) — the form of 
the perturbation H (Si, S2, e) must not be specified in detail. On evaluating ( p6| ) one finds 

- k( ai + faSu + jS 2z ) + \(iik 2 n , (59) 

and u>l ' follows from exchanging the indices (1 «-> 2). The first term is closely related to the classical frequency, 

4 X) S 2z ) = (ft,, ft*) = k dHoi ^ z,S2z) = fc(ai + A5i z + 7 S 2z ) , (60) 

C/Dlz 

and analogously for cjjj 2 ' (Si z , 62*)- In formula (ffi^), the expansion of the generator u p has been given in general 
terms. Focus now on the second term of its right-hand-side containing the difference of the frequency operators in 
the denominator. Its matrix elements in the product basis |mi,TO2}, consisting of eigenstates of the operators Sj Z , 
are given by 

R + R + 
(mi,m 2 |5i + ^^52_|m' 1 ,m 2 ) = (mi, m 2 \-^ L —\m' 1 , m 2 )S mim > i+k d m2m ^ e , (61) 

where, as before, AT, = tD^, 1 -* — £D» 2 \ which in the present case reads explicitly 



(fc(ai + fcS u + JS 2Z ) - £(a 2 + (3 2 S 2z + jS lz )) + ^(fck 2 - f3 2 £ 2 ) . (62) 



Since both operators, Rp k # and A fe ^, depend on the z components only, the expectation value of their ratio, (R/A), is 
equal to the ratio of their expectation value, (R)/(A). Consequently, the convergence of the power series defining the 



first-order generator for the system (57), will be controlled by matrix elements of operators such as in Eq. (J62|) . 
They are made up from two contributions with different origin. The operator in large round brackets is recognized as 
the difference of the frequency operators which one would obtain if in the classical formula ( |60| ) one were to replace 
the variables Sj Z by their operator equivalents, Sj Z . The possibility to approximate arbitrarily well the ratio of 
the frequencies a/ 2 ) / 'a/ 1 ) for a given torus (i.e. fixed values of S\ z and S 2z ) by k/£ leads to the problem of small 
denominators in classical mechanics (cf. [^0[). 

Then, the difference A^ contains a term proportional to Planck's constant: the occurrence of a term linear in h 

follows from studying a Hamiltonian H being quadratic in the z components of the spins. Polynomials of degree r for 
the Hamiltonian lead to terms of at most ?i r_1 , etc. The presence of such terms has been noted already by Robnik 
|PH in a semiclassical treatment of particle systems. He suggested that these additional terms will have the effect to 
move the system away from the classical resonance, however small the value of h. But looking at the exact expression 
( |62| ) for a fixed value of Planck's constant, K = 1, say, one can also conceive it as a new type of resonance condition 
which is no longer linear in the indices k and £: a quantum-mechanical nonlinear resonance condition arises which 
has to be studied in its own right for all possible representations. 

The relevant matrix elements of the operator A^ g (cf. (|6l|)) read explicitly 

(mi — k, m' 2 — £|A^|mi — k, m' 2 — £) = k{a\ + j3\hmi + 7?im 2 ) — £{a 2 + (3 2 hm! 2 + jhmi) — — {flik 2 — (3 2 £ 2 ) , (63) 
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(note the sign change of the last term relative to (p2|)). In each individual representation, there is an important 
difference compared to classical mechanics. The range of the numbers k and I is not all the integers but it varies 
over a finite subset of them only, which is determined by the spin length s. This is characteristic of finite-dimensional 
representations of spin systems while it does not occur for particle systems. The matrix elements of the other operators 
in Eq. (^) lead to expressions of the same structure. 

There are two ways to study the expression j63|) . First, consider it as a quantum- mechanical statement in its own 
right. This means to give H its actual value and to let vary both pairs (mi,m' 2 ) and (k,£) over the finite range 
determined by the value of s. Within each (2s + l)-dimensional representation of the spin algebra, the modulus 
of (Ar«) will have some positive nonzero minimal value, A(s). Consequently, for sufficiently small values of e, the 
expansion of the generator u p in ( |37| ) is expected to converge because the terms (e/A(s)) p can be kept small. For 
spins of length s the perturbed Hamiltonian H e is given by a finite-dimensional hermitian matrix which guarantees 
the existence of a unitary transformation diagonalizing it. Therefore, the generic convergence of the construction 
should not come as a surprise (for possible technical difficulties, cf. fll2|] ). For the present concept of integrability, 
however, it is essential that not one specific representation is studied but that the statements hold in an algebraic 
sense, or, equivalently, in all representations. The implications of this requirement will be made explicit now. 

Imagine to have constructed the diagonalizing transformations in all representations. This set of matrices may 
either be the representations of one single algebraically defined operator U e or not. In the latter case, the unitary 
transformations in the individual representations will not "converge" towards one specific unitary operator for large 
values of s (since the value of h is kept fixed, this procedure is not identical to the classical limit to be discussed 
momentarily). In addition, for larger values of s the number A(s) takes on smaller and smaller values which requires 
the expansion parameter to be correspondigly smaller. If there is no finite value of e such that the perturbation series 
converges for all allowed values of s, i.e. algebraically, the perturbation H (S, e) renders the system non-integrable. A 
strong decrease of the smallest occurring denominator A(s) as a function of s has been observed numerically. If, on 
the other hand, the summation of the perturbation series for all s does not require vanishing e, one may expect an 
algebraic diagonalizing operator to exist; hence, the perturbation can be absorbed into an appropriate redefinition of 
the action operators. In classical mechanics, this convergence is known for finite e is known to result from appropriately 
decreasing numerators. 

Turn now to the second way to study (p3|): in the semiclassical limit with h ~ 1/s and s — > oo, one should recover 
the behaviour of the classical frequency denominators from Eq. (|6^). Numerical tests make it plausible that, indeed, 
the value of A(s) decrease monotonically when larger values of k and t provide better and better approximations of 
the classical frequency ratio (the values of Hmi and Tim! 2 are kept fixed here defining thus a specific "torus" ) . An 
analytic study is not straightforward, even for the quadratically modified resonance condition associated with (|63|): 
the behaviour of the last term is not easily controlled since Tik 2 and M 2 may have values of the order of s. 

VI. DISCUSSION 
A. Summary and outlook 

For the sake of clarity, the main points of the development given in the previous pages are briefly summarized. 
Starting from the algebraic skeleton necessary to describe quantum (spin) systems, arguments for the existence of 
quantum-mechanical action and frequency operators in systems with a single degree of freedom have been put forward. 
When applying these concepts to systems with more degrees of freedom they lend themselves to a definition of quantum 
integrable systems. An algebraic version of perturbation theory of such systems has been formulated capable to 
absorb the effect of an added disturbance. Since this formulation respects (and exploits) the underlying Lie-algebraic 
structure, a tight analogy to classical Lie transforms is achieved. As a natural consequence, operator-valued analogs 
of the classical frequency denominators show up when constructing a unitary transformation which would diagonalize 
the Hamiltonian algebraically. In this way, the quantum-mechanical locus of the convergence problems known to exist 
in classical mechanics has been found, at least on a formal level. In an appropriate sense, the quantum-mechanical 
expressions for the denominators generalize the classical ones. The condition for a classical resonance to occur is, 
effectively, equivalent to a linear relation among the classical frequencies (associated with a "plane" in action space) 
whereas, in a given representation, the relevant quantum-mechanical relation turns out to be nonlinear (defining 
"ellipsoids" or other curved surfaces). Such structures have been observed in semiclassical treatments, for example, 
in the work by Berry and Tabor jl3| who derive energy-level statistics for classically integrable quantum systems in 
the limit of small values for h. 

When scrutinizing the notion of quantum integrability the ultimate goal is to provide evidence for both, its usefulness 
and its consistency. Rigorously testable predictions of specific properties must result for quantum integrable (and 
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nonintegrable) systems, just as the phase-space foliation of classically integrable systems flows from its classical 
counterpart. In our view, the major achievement of the present approach is to have boiled down the question 
of quantum integrability (and, a fortiori, of quantum chaos) to a technical problem, namely to investigate the 
(convergence) properties of a perturbation series. Due to the nonlinearities in the quantum-mechanical expressions 
just mentioned, the behaviour of the perturbation series is, from the present point of view, not obvious. In principle, 
there is no need for quantum systems to exactly parallel the properties of classically (non-) integrable systems: the 
underlying unifying structure is on an algebraic level but with the algebras being realized in different spaces, in 
phase space and in Hilbert space, respectively. In the semiclassical limit, however, one expects that the statements 
associated with the two realizations become comparable if appropriate tools such as Wigner functions are employed. 
The present formulation now calls for a detailed parallel investigation of the effects which a perturbation has on a 
classically integrable system and on its quantum counterpart. 

Here is a brief sketch of a scenario being plausible in view of the present results which, however, has to be substanti- 
ated by further investigations. Focus on the fully quantum-mechanical regime, that is, fix the value of Ti and consider 
the behaviour of the perturbation series in representations of variable dimension, 2s + 1. It is expected that the 
"visibility" of the perturbation will increase for dimensions with larger values of s. A vivid illustration of this effect 
has been reported for a system consisting of two coupled spins where a "regular quantum web" is gradually torn apart 
by a perturbation of increasing strength jL4|. Related phenomena have been observed in particle systems by Bohigas, 
Tomsovic, and Ullmo ]l5j . In a study of the spin-boson model the observed "resonance phenomena" can be attributed 
to the presence of small denominators becoming more or less effective in representations of different dimensions [fl6| . 
An analytic study of these numerically observed effects in the vein of the present paper is under way . In general, 
"quantum resonances" in a nonintegrable quantum system due to ever smaller denominators will become more and 
more pronounced if larger representations of the spin algebra are considered. The infinite-dimensional representation 
(not to be confused with the semiclassical limit) contains the most detailed structure, possibly leading to an actual 
divergence of the perturbation series. 



B. Related work 



There is a large number of joints between earlier work and the present one which so far have been alluded to only 
occasionally. It is useful to dwell on them for a moment, to draw parallels and to point out differences. 

The introduction of frequency operators (in the present terminology) provides a link to the most advanced papers of 
the "old quantum mechanics." In 1925, Dirac Q defined "q-numbers" which would fulfill the fundamental commuta- 
tion relations of position and momentum. Basic algebraic consequences are derived, partly by exploiting the analogy 
to the classical Poisson brackets. One section of the paper is devoted to "multiply periodic systems" characterized 
by the existence of "uniformizing variables" (which nowadays would be called action and angle operators) such that 
the Hamiltonian depends on the action operators only (cf. Eq. @). In the terminology of the present work multiply 
periodic systems are thus recognized as quantum integrable ones. Assuming the action-angle operators to fulfill the 
same commutation relations as momentum and position — which is untenable due to the defectiveness of a phase oper- 
ator jl8) — Dirac introduced two different types of frequencies. The first one is supposed to govern the time evolution 
of the phase operator itself whereas the second one is associated with the time derivative of the exponentiated phase 
operator. It turns out that the latter frequency operator is equivalent to the expression given in (|2p|). Dirac also 
mentions its two possible versions reflecting the freedom of ordering in (|T8|). For spin systems considered here, the 
operator S+ plays a role similar to the exponentiated angle operator. 

Wentzel |19|| continued Dirac's work by establishing a link between Heisenberg's matrix mechanics and the q- 
numbers. In particular, he wrote down expressions for quantum-mechanical frequencies similar to ( |20| ) being applicable 
to multiply periodic systems only. However, this line of reasoning seems to have come to an abrupt end by the advent 
of wave mechanics in 1926: the serious shortcoming of the earlier quantum theory to allow for quantization of multiply 
periodic systems only was removed by Schrodinger's equation which is easily written down for particles subjected to 
arbitrary potentials. 

Similarly, strong formal affinities to Heisenberg's formulation of perturbation theory ]20[ ] exist. This is easily under- 
stood if one recalls that before the introduction of Schrodinger's equation classical mechanics served much more than 
later as Ariadne's thread for quantum-mechanical developments. Still, Heisenberg's approach to treat perturbations 
and the widely used perturbation scheme of Rayleigh-Schrodinger ]2l[ ] differ conceptually in an important way from 
the present one. Traditionally, any quantum system may serve as a starting point for perturbation theory, the only 
assumption being that its (necessarily periodic) solutions are known. In the approach developed here, the unperturbed 
system is assumed to be quantum integrable, hence it presumably stems from a restricted class of systems, and the 
effect of the perturbation is to remove it from this class. In order to prove the KAM-theorem, classical perturbation 
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theory is developed relative to integrable systems: the fate of a torus under the perturbation is investigated, not 
its effect on an arbitrary periodic solution — which might happen to be an isolated periodic orbit of a nonintegrable 
system. KAM-theory would not apply at all. It is hoped that by systematically moving away from integrable quantum 
systems generic properties will be seen to emerge for nonintegrable ones. 

As for a close-up of the relation between classical tori and Heisenberg's matrix mechanics which fits well into the 
scheme developed here, see the work by Greenberg, Klein and Li |2^] . 

Next, there are close ties to Birkhoff normal forms and to the so-called "algebraic quantization" reviewed in 
|p3| . The idea is to quantize integrable approximations to classically nonintegrable systems. Crehan [24f| worked out 
quantum-mechanical normal forms in the spirit of Birkhoff and Gustavson, starting from the classical Hamiltonian 
function in the neighbourhood of a point of stable equilibrium. Impressive agreement between exact quantum- 
mechanical solutions and the approximate ones has been obtained |25|,^6]]. In contrast to these works, however, no 
recourse has been made to classical mechanics in the present approach. 

Concerning the presence of the quantum equivalents of the small denominators, their occurrence in the scmiclassical 
limit has been noted already by Robnik jy] (and by Gram, cf. the introduction of |^|). Robnik pointed out that 
for small values of Planck's constant the classical denominators will be modified by terms of the order T%. It has been 
conjectured that in this way the resonance conditions are destroyed which govern the fate of an individual classical 
torus, that is, its survival or destruction. It is expected that in view of the exact expression ( [f9|) the effect of the 
denominators can be studied in more detail. 

Also, the understanding of clashes between the behaviour of classical and quantum perturbation expansions may 
benefit from the present approach. It has been noted by Eckhardt |28| for a one-dimensional particle system that the 
convergence of a classical perturbation series is not automatically reflected in the convergence of the corresponding 
quantum perturbation series as follows from a simple scaling argument. 

In order to deal with issues of perturbation series for quantum systems a new setting has been provided by Bel- 
lissard and Vittot |^7j whose work is based on the concept of non-commutative geometry. Their work is focussed 
on Lie-transform techniques for quantum-mechanical particle systems with respect to the classical limit. The as- 
pect of integrability, seemingly, has not been addressed in detail. It is particularly important that estimates of the 
Nekhoroshev-type become available which allow one to control the approximations for energy eigenvalues extremely 
well. For systems with one degree of freedom, the convergence of the perturbation series has been proven (cf. (2{|). At 
pre sent, the theory has been worked out for particle systems only while the transfer to spin systems is in preparation 
[po|| . A common feature of many works in this area (cf. also |pl[ is that the unperturbed systems considered are 
collections of harmonic oscillators. Having frequency operators at one's diposal, this restriction seems to be no longer 
necessary. 

Finally, a review of the use of Lie-transform techniques and KAM-like results in quantum mechanics has been given 
by Jauslin [|32|| dealing to a large extent with externally driven systems. 

A: THE COMMUTATORS [F{Sz),S£} 

It will be shown by using the commutation relations of the spin operators that 

[J(S Z ),S\] =Sl(f(S z + kh)-f(S z )) = Sl(exp(khd z )-l)f(S z ), (64) 

where d z is a shorthand for a (formal) derivative with respect to S z . The function f(x) is assumed to have an 
expansion as a power series in x. 
First, it is proved by induction that 

[S* S+] = S+ ( exp(hd z ) - l) S n z . (65) 

This relation holds for n = 1, 

[S z , S+] =S+( exp(hd z ) -ljS z = S + (S z + U-S z ) = hS+ . (66) 

Eq. ( |65| ) is valid for n + 1 if assumed to hold for n: 

[S™ +1 , 5+] = [S z , S+]S Z + S z [S z , 5+] 

= S+ (hS™ + (S z + h)(exp{hd z ) - 

= S+(exp(hd z )-l)S z i+1 , (67) 
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since (S z + K) exp(hd z ) = exp(hd z )S z . Using f(S z ) = J2n=o fn S z implies: 

OO OC 

[/(&), s+\ = J2 f»&> § +] = s+ ( cx p(^) " x ) E ^s n z 

n=0 n=Q 

= S+(exp(hd,)-l)f(S,). (68) 
Thus, Eq. ( fiil ) holds for fc = 1, as is necessary for a proof by induction on k. The step from k to k + 1 reads 
[/(£), = S* [/(&), 5+] + [/(£), 

= ^ +1 ((exp(M z ) - 1) + (exp(fc?i9 2 ) - l)exp(ft0,))/(S,) 

= Sl +1 (exp((k+ l)hd z ) - l)f(S t ), (69) 



where, in the second line, Eq. ( p8|) has been used in the form: f(S z )S + = S + exp(hd z )f(S z ), concluding the proof of 

Relation (|6^) can also be written as 

f(S z )Sl = S h + f(S z + kh), or f(S z -kh)Sl=S k + f(S z ), (70) 

in words, the order of a function of S z only and the A:th power of the operator S + may be exchanged if the argument 
of the function / is shifted by an appropriate multiple of h. Similar relations are obtained for the step-down operator 
by taking the hermitian conjugate of Eqs. (|4|) and (fZOj) . 

A proof of Eq. (68) for arbitrary functions f{S z ) is easily given if one exploits the properties of the ladder operators 
in explicit representations. The action of the step- up operator is in all representations of the form 

S+\m, s) — c s m \m, s) ,c s m = y/ s(s + 1) - m(m + 1) . (71) 

Using f(S z )\m, s) — f(mh)\m, s), one obtains for all |m, s) that 

f(S z )S+\m, s) = C/((m + l)h)\m + 1, s) (72) 
= S+f((m + l)h)\m + 1, s) = S + f(S + h)\m, s) , (73) 

implying that 

f(S z )S+ = S+f(S z + h) , (74) 
which is equivalent to (|68|). As before, induction on k completes the proof of relation (|64|). 



B: CLASSICAL PERTURBATION THEORY 



Results of classical Lie-perturbation theory are collected, for easy comparison with the quantum-mechnical formulae 
an analogous notation is used. Consider a Hamiltonian H (S, e) = Hq(S z ) + H(S, e), where Hq(S z ) is integrable and 
H(S,s) is the perturbation. The goal is to remove the perturbation by introducing a new set of spin variables 

S' = TS, (75) 

where the canonical transformation is implemented by the operator 

? = exp(-{ U (S, £),•}). (76) 

The function u(S,s) is assumed to have an expansion in powers of e, and u(S,e = 0) = 0. When expanding the 
old and new Hamiltonian functions in powers of e, the construction of the transformation T leads to a hierarchy of 
equations 

{H (S z ),u p (S ± ,S z )} = R p (S ± ,S z )-H p (S z ) p = 0,1,2,..., (77) 
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which is the classical counterpart of Q39j) . They can be solved by an appropriate choice of the functions u p and7i p . 
If in classical mechanics the expansion of functions F(S) is written in analogy to (84) the solutions of ( [r7] ) are 



'Hp(S z ) — Rp,o{S z ) , 

,(s±, s.) = £ fsi^4!4 - ^&ls k 



k=l 



iuJk(S z ) iuJk{S z ) 



(78) 
(79) 



The expression for u p (S±, S z ) can be simplified by using tok = from (|l3|). 
Similarly, the results for a two-spin system assume the form of the expansion (|8 



~H P (S Z ) — R P ,oo(S z ) , 

oo / 

u p (S±,S z ) = J2 \ S "+ S ' 



(80) 



fe,i— o 
(fc, 0^(0,0) 



'^((^(S-J+wP^S,)) 



+ l (J k 1 \s z )-J l 1 \s z )) 



S l 2 _ + c.c. 



(81) 



C: EXPANSION OF A FUNCTION F(S±,S Z ) 

A convention for the expansion of functions of spin operators is established. For simplicity, the case of a single spin 
S is considered first. 

Let F(S±,S Z ) be expandable in powers of S+ and S-. Functions depending on S z only are denoted by curly 
symbols J-(S Z ). In general, F will consist of a sum of terms, each of which is of the form 

Sf+S^-fiiSzffi+S^-FziSz) ■ ■ ■ (82) 

The exponents rij ; ± are integers greater or equal to zero. Since the product S+SL = g(S z ) depends only on S z , one 
can write for arbitrary integers n± > that 

J=(|*)_ t ifn+ = n_, (83) 

by exchanging 5(S*z) with S± according to the rules (|70|). Thus, expression (|8^) can be cast into one of these forms 

if now n± are defined as the total number of S± operators in (|8^), n± — Y) . nj±. 

Collecting all terms in the expansion of F(S±,S Z ) with the same number of S+ or 5_, the operator F can be 
written as 

oo 

F(S ± ,S Z ) = F (S Z ) + J2 (S k +F+(S Z ) + F-(S z )Sl) , (84) 
fe=i 

which will be considered as standard form. For the sake of generality, the index k is assumed to run over all integers, 
while in a specific representation the sum in ( |84| ) will run over 2s + 1 values only. If F is a hermitean operator then 
Fq is hermitean and 

F£&)=(Fj;(S z j)\ (85) 

The generalization to two or more degrees of freedom is straightforward but the notation becomes more cumbersome. 
First consider a function F(Si±, S% z , S%±,S2 Z ) depending on the components of two spins Si and S2. As before, the 
function F is assumed to have an expansion in powers of the spin components. Since the components of different 
spins commute, F can be brought to the form: 
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oo 



F(S±±, S±z, S2±, S2z) — F 00 (Si z , S22) + ^ (^Si+S l 2+Fki + (Siz, S2z) 



k,l = D 

(k, ()^(o,o) 



+5{° + F A J (Siz, 5 2z )52_ +S , 2+-F 1 fci + (S'i z , S2 Z )Si_ + F kl (Si z ,S2z)Si_S%_^ . 



(86) 



If F is a hermitean operator this expression simplifies to 



00 



F(S 1 ,S 2 )=Foo(S,)+ J2 (^ + f++(S 2 ) + ^ + F+-(S.)^_+h.c.) , 



(87) 



k,l = 
(fc,i)#(0,0) 



since, generalizing (g|| 




(88) 



Functions depending on the components of more than two spins can be expanded similarly. One gets one term 
depending only on the z components and 2 N terms which contain the explicit dependence on the operators Sj+ and 
Sj-. A unique way of ordering would be, in extension of above, to throughout position the step-up operators Sj+ on 
the left and their adjoints on the right of the middle part which depends on the z components only. 



* This work is based on: Th. Gramespacher (Lie'sche Storungstheorie und die Integrabilitat quantenmechanischer Spinsys- 
teme.) Diplomarbeit (Universitat Basel, 1994). 
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